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INTRODUCTION 


An  important  consideration  in  the  design  of  military  aircraft  is  the 
• • + i -,n  nf  +he  infrared  radiation  emanating  from  the  aircraft  engi  • 

rhe" infrared  s gnature  of  the  engine  can  be  controlled  through  proper - design 
"f  an  engine  diffuser;  hcever,  great  care  must  be  taken  to  assure  that  the 
Proposed  diffuser  does  not  adversely  affect  engine  performance,  iheref  , 
the  iesign  engineer  is  faced  with  the  complex  problem  of  designing  an  e . 
whLb  provides  minimum  turbine  back  pressure  with  an  efficient  exhaust 
diffuser  in  order  to  maximise  shaft  horsepower  and  at  the  same  time  minim, 
radiation  through  the  use  of  curved  wall  diffusers  and  cooled  walls. 

The  satisfaction  of  these  sometimes  conflicting  requirements  P™™" 

to  be  extremely  difficult  in  the  past,  and  engineers  often  have  been  fore 
to  rely  on  empirical  design  methods  based  on  correlations  of  limited 
experimental  data.  For  example,  diffuser  performance  maps  based  on  ‘ 
correlations  have  been  published  by  Reneau  (Ref.  l)  for  ineompressi 
dimensional  flow  and  by  Sovran  (Ref.  2)  for  incompressible  annular  flow  m 
straight-wall  diffusers.  Regions  of  stall  on  these  performance  maps  have  been 
defined  by  Fox  (Ref.  3)  for  two-dimensional  diffusers  and  by  Howard  (Ref.  4) 
for  straight -wall  annular  diffusers.  In  addition,  Sovran  (Ref.  2)  has 
developed1 2 3 4 5 6 empirical  correlations  for  the  effect  of  inlet  blockage  on  perfor- 
mance, and  Runstadler  (Ref.  5)  has  developed  correlations  for  the  effec 
inlet  Mach  number.  Although  these  empirical  design  criteria  provide  some 
insight  into  the  effect  of  variables  such  as  area  ratio,  length,  inle 
blockage  and  Mach  number  on  performance,  these  criteria  are  not  adequate 


(1)  Reneau,  L.  R.,  J.  P.  Johnson,  and  S.  J.  Kline:  Performance  and  Design 

of  Straight,  Two-Dimensional  Diffusers.  Transactions  of  ASME,  Journal  of 
Fluid  Mechanics,  Vol.  89,  March  1969,  PP*  l4l  -l60. 

(2)  Sovran,  G.,  and  E.  D.  Klomp:  Optimum  Geometries  for  Rectilinear 

Diffusers.  Fluid  Mechanics  of  Internal  Flow,  Elsevier  Publishing  Co.,  1967. 

(3)  Fox,  R.  W.,  and  S.  J.  Kline:  Flow  Regime  Data  and  Design  Methods  for 

Curved  Subsonic  Diffusers.  Journal  of  Basic  Engineering,  Transactions  of  the 
ASME,  Series  D,  Vol.  84,  No.  3,  September  1962,  pp.  303 "312. 

(4)  Howard,  J.,  H.  Henseler,  and  A.  Thornton-Trump : Performance  and 

Flow  Regimes  for  Annular  Diffusers.  ASME  Paper  67,  WA/FE-21,  1967- 

(5)  Runstadler,  P.  W.,  and  R.  C.  Dean:  Straight  Channel  Diffuser 

Performance  at  High  Inlet  Mach  Numbers.  Transactions  of  ASME,  Journal  of 
Basic  Engineering,  Vol.  91,  September  1969,  PP*  397-422. 


6 


designing  curved-wall  IR  suppressing  diffusers.  For  curved-wall  annular 
diffusers,  only  a few  studies  such  as  those  by  Dietz  and  Thompson  (Ref.  6) 
and  Thayer  (Ref.  7)  are  available  to  provide  design  information.  In 
particular,  Thayer  has  developed  some  general  design  requirements  for 
diffusers  of  this  type  through  his  investigation  of  the  effects  of  swirl  and 
Mach  number  on  diffuser  performance. 

The  development  of  analytical  design  methods  has  generally  lagged  behind 
empirical  design  methods.  Conventional  solutions,  such  as  those  used  by 
Sovran  (Ref.  2),  divide  the  flow  field  into  an  irrotational  free-stream  flow 
and  a boundary  layer  flow.  These  methods  which  divide  the  flow  field  into 
viscous  and  inviscid  portions  require  an  iteration  between  the  potential  flew 
pressure  field  and  the  boundary  layer  displacement  thickness.  This  iteration 
frequently  fails  to  converge  when  the  boundary  layers  comprise  a significant 
portion  of  the  total  flow  field.  In  addition,  these  iterative  methods  cannot 
account  conveniently  for  phenomena  such  as  inlet  swirl  and  inlet  flow  distor- 
tion. Recently  Anderson  (Refs.  8 and  9)  introduced  a new  method  for  solving 
the  swirling  diffuser  flow  problem  which  solves  a single  set  of  equations 
of  motion  for  the  entire  flow  field  in  the  diffuser,  thereby  enabling 
compatibility  between  the  inviscid  flow  and  boundary  layer  to  be  achieved 
without  the  need  for  matching  a boundary  layer  solution  to  an  inviscid  flow 
solution  through  an  iterative  procedure.  The  method  has  shown  good  agreement 
between  theory  and  experiment  for  incompressible  flow  (Ref.  9)  and  has  been 
extended  more  recently  by  Anderson  (Ref.  10)  to  the  prediction  of  compressible 
flow.  Theoretical  predictions  again  have  been  in  good  agreement  with 
experimental  data. 


(6)  Dietz,  A.  E.,  and  J.  F.  Thompson:  Advanced  Experimental  Infrared 

Energy  Suppression  System  for  the  T-53-L-11  or  T-53-L-13  Turbine  Engine. 

Hayes  Internal  Report  No.  1172,  1968. 

(7)  Thayer,  E.  B.:  Evaluation  of  Curved -Wall  Annular  Diffusers.  ASME 

Paper  71-WA/FE-35}  September  1972. 

(8)  Anderson,  0.  L.:  A Comparison  of  Theory  and  Experiment  for 

Incompressible,  Turbulent,  Swirling  Flows  in  Axisymmetric  Ducts.  AIAA  Paper 
No.  72-42,  10th  Aerospace  Sciences  Meeting,  January  1972. 

(9)  Anderson,  0.  L. : Numerical  Solutions  of  Incompressible  Turbulent 

Swirling  Flows  Through  Axisymmetric  Annular  Ducts.  United  Aircraft  Research 
Laboratories  Report  No.  H2l3577_lj  March  1968. 

(10)  Anderson,  0.  L.:  User's  Manual  for  a Finite -Difference  Calculation 

of  Turbulent  Swirling  Compressible  Flow  in  Axisymmetric  Ducts  With  Struts. 
United  Aircraft  Research  Laboratories  Report  L911211-1,  Contract  No.  NAS3- 
15402,  1972. 
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The  method  derived  in  Ref.  10  requires  construction  of  a generalized 
orthogonal  coordinate  system  from  a solution  of  the  plane  potential  flow 
through  the  duct  in  question.  This  potential  flow  solution  serves  as  an 
approximate  streamline  coordinate  system  upon  which  the  viscous  soluti 
the  equations  of  motion  is  based.  The  equations  of  motion  are  written  in  the 
approximate  streamline  coordinate  system,  and  boundary  layer  approximations 
nay  be  male  in  this  new  system  since  the  potential  flow  streamline 
approximate  the  real  streamlines.  With  the  boundary  ^ “pproximatl^s 
the  equations  of  motion  reduce  to  a set  of  parabolic  partial  differentl  1 
conations  which  apply  to  the  flow  field  under  investigation.  In  Ref.  10,  _ 
r:rton  to  the  potential  flow  problem  was  obtained  through  an  approxi- 
mate geometric  construction  which  yielded  good  results  when  the  curvature  on 
both  walls  was  small  and  nearly  the  same.  This  geometric  solution  sometimes 
failed  when  the  curvature  varied  significantly  from  wall  to  wall.  Even  when 
the  solution  did  not  fail,  significant  errors  could  arise  due  to  the 
curvature  approximation  inherent  in  the  method.  Thus  the  procedure  o Ref.  10 
is  limited  in  the  types  - =i  es  to 

c;“\;;:^taionn:  w:fd::eloPped.  This  new  -tM  obtains  the  solution 
to  (he  plain  potential  flow  problem  using  an  exact  numei  -1  solution  based 
on  the  Schwartz -Christoffel  transformation  (Refs.  11  and  Id). 

It  is  possible  to  solve  the  governing  equations  by  an  explicit  or  an 
implicit  numerical  integration.  In  an  explicit  method  or  so  vt"g  t 
tions  of  motion,  the  allowable  streamwis.  step  size  is  related  to  the  tr 
verse  step  size  through  numerical  stability  conditions.  If  the  solution 
to  b numerically  stable,  a finer  transverse  grid  requires  a smal  er  stream- 
wise  st"  cite.  This  restriction  is  particularly  troubles^  in  the  cas  of 
slot  cooled  walls,  where  a very  fine  transverse  grid  de^ed  *o 
coolant  film  accurately.  Therefore,  under  the  present 
numerical  integration  technique  of  Ref.  10  was  replaced  by  an  implicit 
technique  based  on  the  method  of  Keller  (Ref.  13).  In  this  new  method,  the 


(11)  Kober , H.s  Dictionary  of  Conformal  Representations.  Dover 

Pubiications^r^  Konstruktive  Methoden  der  Konformen  Abbildung, 

Sorineer  Tracts  in  Natwal  Fhilosophy,  Vol.  8,  1963- 
P (13)  Keller,  H.  B.,and  T.  Cebeci  : Accurate  Numerical  Methods  for 

Boundary  Layer  Flows-II  iVo-Dimensional  Turbulent  Flows.  A^  9 
Aerospace  Sciences  Meeting,  New  York,  January  25-27,  1971,  AIAA  Paper 

71-164. 
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equations  of  motion  are  linearized  in  such  a way  that  an  iteration  is  not 
required  to  obtain  a solution.  The  method,  however,  retains  importan+ 
features  such  as  second-order  accuracy  and  from  the  point  of  view  of  linear 
stability  analysis  has  no  restrictions  on  step  size  in  either  the  streamwise 
or  the  transverse  directions  (Ref.  l4).  However,  in  practice  the  step  size 
is  restricted  by  the  need  to  minimize  truncation  errors  arising  from  the 
finite-difference  scheme.  Truncation  errors  will  cause  loss  of  accuracy  and 
may  lead  to  numerical  instabilities  through  nonlinear  effects.  In  addition, 
powerful  matrix  inversion  methods  are  available  in  the  numerical  solution 
(Refs.  15  and  l6). 

Under  the  present  effort,  advanced  mathematical  techniques  have  been 
applied  to  the  solution  of  turbulent  compressible  swirling  flow  through 
curved -wall  annular  diffusers  with  slot-cooled  walls.  From  this  analysis,  a 
computer  program  has  been  developed  and  sample  cases  calculated  and  compared 
with  experimental  test  results  obtained  from  an  IR  suppressing  diffuser  with 
slot -cooled  walls  at  different  simulated  engine  operating  conditions. 


(14)  Keller,  H.  B. : A New  Difference  Scheme  for  Parabolic  Problems. 

Numerical  Solution  of  Partial-Differential  Equation -II  SYNSPADE  1970  Ed.  by 
Hubbard,  B.  Academic  Press,  New  York. 

(15)  Keller,  H.  B.:  Accurate  Difference  Methods  for  Linear  Ordinary 

Differential  Systems  Subject  to  Linear  Constraints.  SIAM  J.  Namer,  Anal. 
Vol.  6,  No.  1,  March  1969. 

(16)  Briley,  W.  R.,  and  H.  McDonald:  An  Implicit  Numerical  Method  for 

the  Multidimensional  Compressible  Navier-Stokes  Equations.  United  Aircraft 
Research  Laboratories  Report  M91l363-6,  November  1973- 
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ANALYSIS 


The  present  analysis  solves  the  problem  of  axisymmetric  swirling  flow 
through  typical  IR  suppressing  diffuser  geometries  in  a two-step  procedure. 

In  the  first  step,  a proper  coordinate  system  is  constructed;  in  the  second 
step,  a set  of  boundary  layer  type  parabolic  partial  differential  equations 
is  solved  using  a forward  marching  implicit  numerical  integration  procedure. 

For  flow  over  a flat  plate,  the  proper  coordinate  system  consists  of  lines 
parallel  to  the  plate  (termed  the  streamwise  coordinate)  and  a second  set  of 
lines  perpendicular  to  the  plate  (termed  the  transverse  or  normal  coordinate). 
If  the  equations  of  motion  are  written  in  this  Cartesian  coordinate  system  and 
the  boundary  layer  approximations  are  made,  a set  of  parabolic  partial 
differential  equations  is  obtained.  For  this  simple  problem,  it  is  obvious 
that  the  boundary  layer  approximations  (namely,  that  the  transverse  velocity 
is  small  compared  to  the  streamwise  velocity  and  that  the  streamwise 
derivatives  are  small  compared  to  the  transverse  derivatives)  are  valid.  In 
the  case  of  more  complicated  geometries  such  as  curved-^wall  diffusers,  the 
coordinate  system  in  which  the  boundary  layer  approximations  can  be  made  is 
not  as  simple  as  in  the  Cartesian  coordinates  described  above.  Rather,  it  is 
a coordinate  system  in  which  one  coordinate  approximates  the  streamlines  and 
the  other  is  normal  to  the  streamlines.  Such  suitable  coordinates  can  be 
obtained  from  the  plain  potential  flow  solution  for  the  duct  under 
investigation  since  it  is  apparent  that  in  view  of  the  constraining  effect 
of  the  walls,  the  potential  flow  streamlines  approximate  the  real  streamlines 
provided  large  regions  of  flow  separation  do  not  occur. 

This  problem  has  been  discussed  in  more  detail  by  Anderson  (Ref.  10), 
where  it  was  shown  that  the  solution  to  the  plane  potential  flow  problem 
through  a given  duct  can  be  used  to  construct  an  orthogonal  coordinate  system 
uniquely  suited  to  solve  for  the  turbulent  flow  through  the  duct.  Although 
the  direct  problem  of  determining  the  velocity  potential  s and  stream  function 
n in  terms  of  the  cartesian  coordinates  R and  Z may  be  solved  easily,  the 
equations  of  motion  for  the  turbulent  flow  require  that  n and  s be  explicitly 
the  independent  variables  so  that  the  coordinate  functions  R(n,s)  and  Z(n,s) 
can  be  obtained  (Ref.  10).  Although  an  approximate  solution  to  this  problem 
was  presented  in  Ref.  10,  the  solution  is  inaccurate  for  ducts  having  large 
curvature.  In  order  to  alleviate  this  curvature  limitation,  an  exact 
numerical  solution  has  been  obtained  to  the  inverse  problem  using  the 
Schwartz -Christoffel  transformation  (Ref.  11).  In  addition,  the  Schwartz- 
Christoffel  transformation  provides  a method  for  obtaining  an  orthogonal 
coordinate  system  for  a duct  with  slots.  The  method  for  solving  for  the 
generalized  orthogonal  coordinate  system  is  derived  in  the  next  section. 
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The  boundary  layer  approximations  to  the  equations  of  motion  for 
turbulent  swirling  flow  with  normal  pressure  gradients  are  derived  in  Ref.  10 
where  these  equations  were  solved  using  an  explicit  numerical  integration 
method.  For  slot -cooling  problems,  however,  this  explicit  method  is 
unsuitable  because  the  inner  layer  of  the  turbulent  boundary  layers  cannot  be 
described  accurately  if  a reasonable  streamwise  step  is  to  be  taken.  There- 
fore, under  the  present  effort,  an  implicit  method  of  numerically  integrating 
the  equations  of  motion  was  developed.  This  method  is  derived  in  the  follow- 
ing section  entitled,  "Conformal  Mapping  Solution".  The  implicit  method  is 
unconditionally  stable,  and  the  linearization  technique  used  permits  integra- 
tion of  the  equations  of  motion  without  any  iteration  such  as  that  used  by 
Keller  (Ref.  13). 

Conformal  Mapping  Solution 
Schwartz -Christoff el  Transformation 

If  a curved-wall  duct  is  represented  by  straight  line  segments  in  the  w 
complex  plane  to  form  a many  sided  polygon,  the  Schwartz -Christoffel  trans- 
formation (Ref.  11)  may  be  used  to  transform  this  polygon  in  the  w plane 
into  the  upper  half  of  the  z plane,  as  shown  in  Fig.  1.  Under  this  trans- 
formation, the  source  flow  at  the  duct  inlet  in  the  physical  plane  becomes 
a point  source  at  the  origin  of  the  z plane.  Source  flows  resulting  from 
inlet  cooling  slots  become  point  sources  on  the  real  axis  of  the  z plane. 

The  potential  flow  solution  as  a result  of  this  source  distribution  in  the  z 
plane  can  be  found  easily  by  superposition  of  elementary  source  solutions 
leading  to  a definition  of  the  streamlines  n and  potential  lines  s in  the  z 
plane.  Then  given  n and  s in  the  z plane,  R(n,s)  and  Z(n,s)  can  be  obtained 
by  going  back  to  the  w plane  and  rotating  and  scaling  as  shown  in  Fig.  2. 

This  procedure  is  explained  in  more  detail  in  the  following  paragraphs. 

The  conformal  mapping  method  has  several  important  advantages  over 
other  methods  of  determining  R and  Z as  a function  of  s and  n.  First,  the 
inverse  problem  can  be  solved  exactly  in  a straightforward  manner  as  opposed 
tc  most  other  procedures,  which  lead  to  approximate  solutions.  Second,  real 
ducts  do  have  discontinuities  along  the  wall  boundaries  for  which  the 
Schwartz -Christoffel  transformation  is  ideally  suited.  Third,  the  technique 
developed  in  this  report  determines  the  wall  slope  and  integrates  the  slope 
to  obtain  the  wall  contour.  Thus  the  first  derivatives  and  metric  scale 
coefficients  required  for  integration  of  the  viscous  flow  equations  are 
obtained  directly  rather  then  by  numerical  differentiation,  leading  to  a more 
accurate  solution. 
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Solution  in  z Plane 


The  complex  potential  for  a source  located  at  the  origin  of  the  z plane, 
which  represents  inflow  at  the  duct  entrance  plane,  is  given  by 


F = l nz  = S + in 


(1) 


The  complex  potential  can  be  solved  explicitly.  Thus,  as  shown  in  Fig.  1, 


s = | nr  = In  ^x2  ♦ y2j 


(2) 


n = <fi  - tan 


(3) 


Equations  (2)  and  (3)  describe  the  potential  flow  in  the  z plane  in  the 
absence  of  any  slot  cooling.  Since  the  walls  contain  poles  representing 
corners  of  the  duct  in  the  physical  plane,  a finite  upper  half  plane  bounded 
by 


( <<fi  S w -£ 


(4) 


(5) 


is  defined.  Thus  n and  s are  bounded  by 


£ < n ^ v -£ 


(6) 
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(7) 


In  rQ  £ s < In  rL 


and  the  solution  lies  completely  within  a bounded  domain  free  of 
singularities . 

Solution  in  w Plane 

4 

The  Schwartz -Christoff el  transformation  is  given  by 

a*  . 4.  r /2-bI'r0l/’'  (8) 

d2  2 1=2  \ 1/ 


where  b-j-  is  the  location  of  the  poles  on  the  x axis  of  the  z plane, 
representing  comers  in  the  physical  plane,  and  aj  is  the  corresponding  corner 
angle  (defined  in  Fig.  l)  in  the  w plane.  The  bj's  and  Qfj's  are,  therefore, 
real  constants.  When  the  values  of  bj  are  known,  any  point  (n,s)  in  the  z 
plane  corresponds  uniquely  with  a point  in  the  w plane.  The  central  problem 
is  to  find  the  values  for  the  bj's  which  are  unique  for  the  duct  under 
consideration.  Let 


W = £ + iz 


(9) 


z = x +iy 


(10) 


Then,  because  of  orthogonality, 


dw 

dz 


.ii-i.4L.fUi.gZL 

dx  <Jy  <Jy  dx 


(11) 
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The  real  and  imaginary  parts  of  Eq.  (8)  are  evaluated  as  follows : 

rj=|(x-bjf  + y‘F/2  m 

^ = ton'1  |y/(x-b,)|  (13) 

Vr/<W’r  W 


v/hich  can  be  evaluated  by  repeated  application  of  the  product  rule  for  complex 
numbers , 


Xj*l  * Kj  xj*i  ~ 


(19) 


7 = x y + y • x 

j*i  j j ♦ i j j* 


(20) 


Finally,  comparing  Eq.  (l8)  with  Eq.  (ll)  results  in 


7 - il  - iZL 
x ' ax  ~ dy 


(21) 


~ _ M. . iH 

^ dy  dx 


(22) 


which  lead  to  differential  equations  relating  §,  T|,  and  x,y.  As  shown 
subsequently,  these  relations  allow  the  construction  of  the  required  potential 
solution  in  the  physical,  w plane. 

Differential  Equations 

The  next  step  in  the  solution  requires  the  derivation  of  the  differential 
equations  valid  along  an  n or  s coordinate.  From  Eqs.  (2)  and  (3), 


as  _ dn_  _ s 

dx  dy  x2  + y2 


(23) 


Al In.  = y_  (24) 

dy  d y x2  + y2 
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A determinant,  D,  is  defined  by 


0 = ' [(-aff+  (W)  1 


(25) 


Then 


■+  l 


- dn-  4t  ds 

a*  ay 


(26) 


dy  = -fr  ["  If  dn‘  "af ds]  = i [' 


an  . , an  . 
ay" dn+  ax  ds 


(27) 


Hence  along  a streamline  dn  = 0 and 


1>L 
a s 


_ I -&S_ 
If  dx 


(28) 


d± 

as 


J_  jL§_ 
D ay 


(29) 


Along  a potential  line  a r ~ 0 and 


ax 

an 


j_  _an_ 
d ax 


(30) 


ay  i an 

an  o ay 


(31) 
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Equations  (29)  through  (31 ) allow  construction  of  the  solution  in  the  z plane. 
Finally,  using  Eqs.  (ll),  (2l),  and  (22),  an  integration  may  be  carried  out 
along  streamlines  or  potential  lines  to  construct  the  solution  in  the  w plane. 


dC  _ dx_  il  . djl 

ds  ’ dx  ds  dy  ds  ’ dn 


(32) 


ds  dx  ds  dy  ds  " dn 


(33) 


Hence,  integration  along  streamlines  in  the  w plane  is  obtained  through  Eqs. 
(32)  and  (33)  together  with  Eqs.  (27)  and  (29),  and  integration  along 
potential  lines  in  the  w plane  through  Eqs.  (32)  and  (33)  together  with  Eqs. 
(30)  and  (31).  The  metric  scale  coefficients  are  the  same  in  both  directions 
and  are  given  by 


v 


d F /dw 
dz/  dz 


(34) 


which  is  the  same  as  the  magnitude  of  the  potential  flow  velocity  obtained 
from  the  complex  conjugate.  It  should  be  noted  that  the  solution  is  equiva 
lent  to  a solution  for  the  wall  slopes,  Eqs.  (32)  and  (33)?  or  the  metric 
scale  coefficient.  Hence  the  wall  slopes  and  metric  scale  coefficients  are 
solved  for  directly  and  the  wall  contour  is  obtained  by  integration.  The 
determination  of  the  bj's  which  define  the  duct  is  discussed  subsequently. 

Solution  Near  a Source 


Consider  the  solution  in  the  neighborhood  of  the  inlet  source  as  z » 0. 
From  Eq.  (8) , 


dw 

dz 


(35) 


where  C]_  is  a complex  constant  given  by 
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(36) 


■ cn  (s-so)  + cm  (n'n°) 


„ and  c are  the  real  and  Imaginary  parts  of  Cr  Iherefore  the 
where  and  are  axis  of  symmetry  (see 

inlet  is  a straight  duct  with  an  angle  \alo 

Fig.  1) 


a0  ■ ton 


'(cIr/Cir) 


From  Eq.  (34),  the  metric  scale  coefficient  is  given  by 

I 2 „ 21  1/2 

V |C|I  + CIB  | 


However,  at  the  inlet,  From  Eq.  (2)  and  Eq.  (3) 


s-s0  = 


(42) 


n-n0  = <*> 


(43) 


and  from  Eq.  (4), 


An  = n - 2£ 


(44) 


Therefore,  inlet  height  is  given  by 


h = AH  (45) 

o v0 


Transformation  to  R,Z  Plane 

Since  the  inlet  of  the  duct  starts  out  with  an  angle  a0,  as  shown  in 
Fig.  1,  the  transform  to  the  (r,z)  plane  shown  in  Fig.  2 is  obtained  through 
a rotation  of  an  angle  aQ  and  is  given  by 

r-r0  = cos  a0  (17-17J  - sina0  ($-£0)  (46) 


2-20=cosa0($-C0)  + ^a0(v-%) 


(47) 
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Finally,  the  transformation  to  the  (R,Z)  plane  is  obtained  through  a 
translation  and  scaling  using  the  inlet  height  h0: 

+ (48) 

Z-.  (Rtq-  V_  W) 

The  streamline  coordinates  are  scaled,  noting  that 

£ < n < n -i 


0 < n < i 


(51) 


0 < S < sL 


(52) 


Thus 


n-C  n-C 
= n-2C  " 


(53) 


s~so  . S So  (54) 

• n-2C  " 
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(55) 


v = 


! _y_ 

rto”rho  vo 


For  numerical  convenience 


V - (t  - 2 0 A/2 


(56) 


Then  from  Eq.  (42) 


ro  = exp  S0 


(57) 


and  the  solution  in  the  z plane  is  located. 

Locating  Poles 

The  location  of  the  poles  bj  in  the  z plane  must  be  obtained  by  an 
iterative  method.  It  is  noted  that  the  location  of  the  corners  and  the 
corresponding  angle  change  of  the  polygon  representing  the  physical  duct  in 
the  (R,Z)  plane  is  known.  Specifically,  the  location  of  these  corners  can  be 
expressed  in  terms  of  the  distance  X(j)  along  the  wall  from  the  inlet  to  the 
jth  corner.  Each  corner  represents  a pole  in  the  z plane.  If  a guess  is  made 
for  the  bj's,  then  Eqs . (32)  and  (33)  “ay  be  integrated  along  the  wall  stream- 
lines. Since  the  location  of  each  pole  is  known  in  the  (R,Z)  plane,  the 
distance  X(J)  to  the  j44  pole  is  computed  and  compared  to  the  known  X(j). 

An  iteration  procedure  based  on  Newton's  method  is  used  to  obtain  a new  guess 
for  the  bj's.  In  the  iterative  procedure  the  duct  contour  is  defined  by 
specifying  the  wall  radius  at  JL  equally  spaced  mesh  points.  Define 


AZ  = ZL/(JL-  I) 


(58) 
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(59) 


Zj  = AZ  (J-i) 


Then  the  hub  and  tip  contours  (outer  wall  and  inner  wall)  in  the  physical 
plane  are  known  at  each  of  the  J points. 


RH(J)  = R 


(60) 


Rt{j)  = rt(zu) 


(61) 


0H  (j)  = tan-1 


(62) 


0T(J)  = tan'1  (-^ ) 


(63) 


The  Oj’s  for  the  Schwartz -Christoff el  transformation  are  then  given  by 


aH  (j)  = eH  0)  - eH  (j  - 1) 


(64) 


aT(j)  = 0T  (j+i)-$t(j) 


(65) 
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where  they  ha/e  been  defined  as  the  change  in  wall  angle  moving  around  the 
duct  contour  (polygon)  in  a counterclockwise  direction.  Since  the  polygon 
is  composed  of  straight  line  segments,  the  distance  along  the  wall  from  the 
inlet  to  the  Jth  point  is  given  by 

[RH(l)-RH(I-l)]2+  Az2  j 1/2  (66) 


XT(J)  = t2  j [rtU>  - r7(I-0]2  + A z2  | 1/2  (67) 


If  an  initial  guess  is  used  for  the  solution  using  the  approximate 
solution  described  in  Ref.  10,  then  the  poles  may  be  located  using  Eq.  (42 ); 
Eq.  (32)  and  Eq.  (33)  may  be  integrated  along  the  walls  (streamlines)  and  the 
distance  along  the  wall  to  each  corner  determined  from 


X*(J) 


(68) 


x> 


(69) 


These  X^(J)  and  X^(J)  in  general  will  not  agree  with  that  calculated  using 
Eqs.  (66)  and  (67).  However,  at  each  pole  in  the  (R,Z)  plane  we  may  obtain  a 
new  guess  for  b-[  using  Newton's  method 

S*  (J)  = SH  (J)  + v>  |x*  (J)  - XH(J)  j ( 70 ) 
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(71) 


x*(j)  - XT(J) 

Convergence  occurs  when 


u*\  V v 

ST(J)  = st(j)  + vT  (J) 


| xtoi- X(J)|  £ {M 


(72) 


v+1 

for  all  J.  Once  the  S(J)  are  known,  the  new  location  of  the  poles  may  be 
obtained  using  Eqs . (42)  and  (54). 

Duct  With  Slots 


Inlet  slots  in  a duct  must  satisfy  the  Kutta  condition  that  the  stream- 
line leaves  tangent  to  the  slot  lip.  Thus  the  Kutta  condition  is  equivalent 
to  stating  that  the  static  pressure  on  each  side  of  the  slot  lip  is  the  same. 
Because  of  the  Kutta  condition,  the  coordinates  for  ducts  with  slots  may  be 
calculated  by  overlaying  solutions  of  successively  larger  ducts  without  slots. 
This  procedure  is  shown  schematically  nn  Fig.  3*  The  coordinates  are 
calculated  for  duct  1 from  station  (l)  to  station  (2).  Then  the  cordinates 
are  calculated  for  duct  2 from  station  (2)  to  station  (3)*  T1 us  the  Kutta 
condition  is  satisfied  at  the  slot  lip  by  construction  of  the  streamline 
(wall)  for  duct  1 tangent  to  the  slot  inlet.  This  process  may  be  repeated 
for  any  .lumber  of  slots  . 


Implicit  Method  of  Solution 
Mach  Number  Transformation 

At  low  Mach  numbers,  the  extremely  small  variation  of  the  pressure  and 
temperature  within  the  diffuser  leads  to  large  numerical  errors  in  the  solu- 
tion of  the  equations  if  the  actual  pressure  and  temperature  are  treated  as 
dependent  variables.  Therefore,  a Mach  number  transformation  was  devised  in 
which  the  dependent  variables  are  the  difference  of  the  local  pressure  and 
temperature  from  the  mean  inlet  flow  conditions.  For  the  purpose  of  the 
transformation,  ll,  0,  I are  definedjis  the  mean  inlet  pressure,  tempera- 
ture, and  entropy,  respectively,  and  II,  9,  I , Q are  defined  by  the 
relations 
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n =n,  + yMr2  FI 


(73) 


®:  ®,  + (y-l)  Mr2  ® 


(74) 


l = V(r-')Mr2T 


(75) 


o = (y-i)  Mr2  5 

where 

r y- 1 In  in  IT, 


(76) 


(77) 


The  variables 
noted  that  for 


__  ~ ~ 

11 » I>  Q are  the  new  dependent  variables, 

very  small  Mach  numbers,  Eq.  (73)  becomes 


It  should  be 


n-n, 

77z  = 0 (')  (78) 

y M 


Whe°  this  transformation  is  applied  to  the  equations  of  motion  as  given  in 
Ref.  10,  all  the  terms  in  the  equations  become  the  same  order,  allowing  an 

te°rT4e  tn"6  S°1Utl0n  t0  te  °btalned-  In  the  P™*'™  Emulation,  some 
terms  in  the  governing  equations  were  eonsiderably  larger  than  others,  leading 
to  numerical  errors.  ’ J-cauj-nK 
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Baric  Equations  of  Motion 

Unier  the  present  effort  the  Mach  number  transformation , Eqs . ( 73) 

' is  applied  to  the  equations  of  motion  derived  by  Anderson  . 

(He^lO  in  addition,  the  eqaations  are  arranged  as  first-order  equations 
lo  facilUate  the  application  of  an  illicit  na»ericai  integration  method. 
The  governing  equations  are: 

Continuity  Equations 


d-n 


=o 


Streanwlse  Stress  Component 


Tangential  Stress  Component^ 


Normal  Momentum  Equations 


an 

dr) 


+ 


_L  4^-lPU 
xr  au,ru 


}4> 


Entropy  Equation 

injlr-OM'S]  -in  [yM*n] 


(79) 


(80) 


(81) 


(82) 


(83) 
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Heat  Flux  Equation 


/ ' /*E\ 

a®  nr 

\PRE  / 

dr)  XV 

(84) 


Equation  of  State 


n,  + yMr2  fl  = Pj®,  + (y-l)  Mr2  £>] 


(85) 


Streamwise  Momentum  Equation 

al 


ns 


V 

IT 


dr) 

d±  _aus 

dr)  dS 


v a / g 

r)-  _L 

. JV. 

xg  an 

r)  XV 

dn 

+ [xl  XL 

aus 

f [JL 

[ g J as 

dr) 

[ XR 

1 an 

f Hs 

x TT  : 

XV 

ns 


dR 


PU 


* 


(86) 


Tangential  Momentum  Equation 

J V ±(1WJ-  ii]y 
drj  [1<^  an  \Wtxr  anj^n^. 

- [XI  il  + [xl  M ^ 
l G I dr)  as  + m as  dr) 

ft] 
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Energy  Jlquati  PH 


(88) 


It  is  noted  upon  examination  of  Eqs.  (79)  through  (88)  that  no  derivatives 
of  density  appear  explicitly;  therefore,  a Mach  number  transformation  was  not 

applied  to  the  density. 

Boundary  Conditions 

For  annular  flow  fields  within  ducts  having  both  an  inner  and  an  outer 
wall,  the  proper  set  of  boundary  conditions  is  (see  Ref.  10): 

^(S,0)  = %(S) 
us(sp)  = o. 

4f>(s,o)  = 0. 

Q (s,o)  = 0.  adiabatic  wall 

or  ® (S.O)  = ®H 


*(S,i)  = *T(S) 

Us(S,l)  = o. 

iys.1)  = 0. 

q (S,D  = o.  adiabatic 
or  ® (S.O  ■ 0T(S) 
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where  S is  the  streamwise  coordinate  and  Tl  is  the  transverse  coordinate. 

The  transverse  grid  is  normalized  so  that  the  walls  occur  at  71=  0 and  7X  = 1. 

For  axi symmetric  flow,  in  which  no  inner  wall  is  present,  the  equations 
of  motion  contain  a removable  singularity  at  the  origin  or  axis  of  symmetry. 
The  boundary  conditions  (Eq.  (89))  must  be  replaced  by  boundary  conditions 
based  on  a Taylor  series  expansion  of  the  flow  variables  about  the  centerline. 
From  Eq.  (86)  at  a small  distance  h from  the  centerline,  the  expansion  for 
LnS  is  given  by 


l 


+ V 


#]0h  + 4s) 


(91) 


The  remainder  is  neglected  because  it  is  of  higher  order  than  the  order  of 
the  difference  approximation.  Then  Eq.  (80)  and  Eq.  (91)  yield 


U.  = u 


so 


NR 


(92) 


which  serves  as  a boundary  condition  for  Us . 

The  expansion  of  En^  about  the  centerline  is  obtained  from  Eq.  (87): 

*V4’)  (93) 

Hence,  from  Eq.  (8l), 

u<*>,  =°(h6)=0  (94) 

since  it  is  of  higher  order  than  the  difference  approximation.  The  heat  flux 
equation  (Eq.  (84))  and  the  energy  equation  (Eq.  (88))  are  used  to  find  a 
boundary  condition  for  Q: 
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(95) 


Q = - 


1 

2 


r-i 

y 


VPUg® 

. 


® = ®0  + T 


r-i 


IvPu.® 

XV  [ * 


(96) 


And,  finally 


, a boundary  condition  for  (//  is  obtained  from  Eq.  (79) 


* = 2ll  (Pus)0  h2 


(97) 


The  boundary  conditions  may  then  be  applied 
line.  However,  a great  simplification  may 


at  a distance  h from  the  center - 
be  obtained  if  h is  chosen  such 


that 


h/  (vAt  j)  « I 


(98) 


which  implies  that  the  first  point  is  very  near  the  centerline.  Then  we  have 
the  axi symmetric  boundary  conditions 

V (s,o)  =o. 


I Js,o)  =0. 

ns 

(S,0)  =0. 

3 ts.o)  =o. 


(99) 
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Finite  Difference  Approximation 


These  equations  are  reduced  to  finite  difference  equations  using  the 
second-order  finite  difference  scheme  of  Keller  (Refs.  14  and  15 ).  The 
following  notation  for  S and  T)  is  introduced: 


ASJ  = SJ  - SJ"' 

Ar?K= 

sJ-|/2  =4,(sJ+  sJ-') 
"Vl/2  = T(\+  ^K-l) 


and  for  any  dependent  variable  g (z,s)  g 2 and  gK  ^ are  defined  by 


(100) 


J-l/2  I ^ 

g - t lQ  +g  ) ( (ioi 

-i/2  = T (V  9k-,)  j 

The  equations  of  motion,  Eqs . (79)  through  (88),  are  linearized  by  performing 
a Taylor  series  expansion  in  the  S coordinate,  as  suggested  by  Briley  and 
McDonald  (Ref.  l6).  Let  any  dependent  variable  g be  given  by 


gJ  = gJ_l+  Ag 


(102) 


where 


(103) 


Then  we  have  the  following  product  rules  (Ref.  16): 


31 


Hj  j j-i  ,j-i  J /#  \J-i 

i = f q ♦ f g - (fg) 

/ \J- 1/2  | J J-I  , J-*  J\ 

(fg)  = i(fg  +f  O 

(fgh)J  = (fg),',hJ  + (gh)J',fJ+(fh),‘lgJ  - 2 (fgh)1 

Ug\J-,/2_  < ~ <1 

lasjK  ‘ as 

/_ag_ ‘ 1/2  gJ-gJ~'  fJ-' 

l ac  f ' TVs-  fK 


Substitution  of  Eq.  (104)  into  the  equations  of  motion,  Eqs . (79) 
through  (88),  yields  the  following  results. 
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Normal  Momentum  Equation 
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Entropy  Equation 
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Heat  Flux  Equation 

(j-  J1\J"  ( 

\Per  fJ'r  ' K —1/2 

§J-sJ 

, K K- 

\ . A7? 
i)+  T 

MJSk+5kh)=0  (110) 

Equation  of  State 

m2  ni  - 
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Streamwise  Momentum  Equation 


3M 


(112) 


Tangential  Momentum  Equation 
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Energy  Equation 
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With  mass  flow  bleed  at  the  wall,  ^(s)  and  s)  are  given  by 


L. 


{4-r*  <«) 

tTJ  =^j"'  ♦ As,  (^L)^  ^ (116) 


Solution  of  Matrix  Equation 


The  solution  of  these  equations  is  obtained  using  the  method  of  block- 
tridiagonal  factorization  (Refs.  15  and  17).  If  the  column  matrix  JK  is 
defined  by 


(117) 


then  the  difference  equations,  Eqs.  (105)  through  (ilU),  may  be  written  as  a 
matrix  equation 


SK  (118) 


where  and  are  the  coefficients  of  the  dependent  variables  and  is  the 
right-hand  side  of  these  equations.  If  T]j_  =0  is  at  the  first  mesh  point  and 
%l  = 1 is  at  the  last  mesh  point,  equations  are  written  at  each  of  the  T)k 
transverse  locations  where  2£K^KL  and  ten  boundary  conditions  are  required. 
Note  that  Eqs.  (109)  and  (ill)  involve  only  the  mes^  point.  Tnen  we  may 
use  these  two  equations  as  boundary  conditions  plus  the  eight  conditions 
given  by  Eqs.  (89)  and  (90)  or  Eqs.  (89)  and  (99)- 


(17)  Varah,  J.  M.:  On  the  Solution  of  Block-Tridiagonal  Systems  Arising 

from  Certain  Finite -Difference  Equations.  Mathematics  of  Computation, 

Vol.  26,  No.  1,  March  1969. 

36 


...... 


A 5X10  M matrix  for  the  boundary  conditions  may  then  be  written  for  each 
end  point  such  that  the  matrix  equations  become 

a'T'=f'  | (119) 

aKLTKL=rKL  1 

The  complete  set  of  matrix  equations  is  written 

A 7 = Q <12°) 


where 


3 ={?',  T\Tkl}t  (122) 
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The  matrix  A la  made  block  tridiagonal  by  splitting  as  follows  : 


(123) 


(124) 


(125) 


(126) 


(127) 
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Hence  we  have  a block  tridiagonal  matrix: 


A1  C ' 

•*2  ¥2  5?  2 


r v e 


BKr£K 


BKL-'  S KL_1 5KL_I 


bkl  zkl 


The  matrix  Q splits  as  follows: 


q'=  ?,*(!•  1, 5)|T 


f=  {t/ (1=6,10),*“  (I- 'JBt} 


5^{TIKLa.6,io),f^} 


The  matrix  equations  (Eq.  (120))  are  solved  using  the  method  of  block  tri- 
diagonal factorization  (Ref.  17)  with  the  recursion  formulas  given  by 


*'=  A1 


I*.  (CV';K 


(133) 


DK=r-iKr‘'  2-k4kl 


(134) 


z‘=  (dT'o1 


z"  (oV  qk-bkzk'' 


2 i K-  KL 


and  the  solution  is  obtained  by  backward  substitution. 


f K=L  zKL 


= K |K_|K^K+I  ( Kl_  — I > Z,-0 


4o 


(135) 
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DESCRIPTION  OF  TEST  FROGRAM 


The  3T9  demonstrator  IR  suppressing  exhaust  diffuser  was  tested  on  the 
Pratt  & Whitney  Florida  Research  and  Development  Center  (FRDC)  D-32  test 
facility  to  verify  the  accuracy  of  the  analysis  described  previously  in  this 
report.  The  FRDC  facility  can  provide  hot  gas  flow  at  a temperature  of  1200 
deg  F and  weight  flows  of  8.3  lb/sec.  Adjustable  swirl  vanes  upstream  of  the 
diffuser  provide  swirling  flow  from  0 deg  up  to  30  deg,  thus  simulating 
typical  turbine  exit  flows  at  different  engine  operating  conditions.  In 
addition,  the  facility  was  modified  to  provide  coolant  flow  rates  up  to  10 
percent  of  diffuser .airflow  rates  and  to  provide  slot  cooling  on  the  diffuser 
walls.  Complete  wall  static  pressure  and  temperature  instrumentation  was 
provided  to  measure  the  effect  of  wall  cooling  rates  on  pressure  recovery  and 
wall  temperature  at  different  simulated  engine  operating  conditions. 

Description  of  Test  Facility 

A schematic  of  the  FRDC  D-32  test  facility  is  shown  in  Fig.  4.  The 
diffuser  inlet  flow  simulating  a typical  turbine  exit  flow  is  provided  by 
compressed  air  from  the  D-Area  J57  slave  engine.  The  compressor  air  flows 
through  a stand  heater  burner  and  swirl  generator  with  adjustable  vanes 
shown  in  Fig.  5.  Coolant  airflow  is  supplied  to  the  inner  and  outer  walls  of 
the  exhaust  diffuser  from  the  D-Area  350-psig  air  system.  As  shown  in  Fig.  6, 
the  coolant  to  the  outer  wall,  inner  wall,  and  base  louvers  is  separately 
regulated.  A schematic  of  the  ST9  demonstrator  IR  suppressing  exhaust 
diffuser  is  shown  In  Fig.  7.  The  manifold  providing  the  coolant  flow 
separately  to  the  inner  and  outer  walls  and  base  is  shown  as  well  as  the 
location  of  the  louvers  or  slots  on  the  diffuser  wall.  Throughout  the  test 
program,  no  coolant  flow  was  supplied  to  the  base  in  order  to  eliminate 
effects  of  base  cooling  flow  on  inner  wall  static  pressure.  Photographs  of 
the  test  facility  with  the  diffuser  in  place  taken  after  completion  of  testing 
are  shown  in  Figs.  8 and  9*  No  apparent  damage  or  deterioration  to  the 
suppressor  occurred  during  the  test  program. 


Description  of  Instrumentation 


Wall  pressures  were  measured  by  thirty-one  static  pressure  taps  located 
on  the  inner  and  outer  walls  of  the  diffuser  with  the  data  recorded  manually 
on  manometer  boards.  The  exact  axial  locations  of  these  pressure  taps, 
denoted  IW01  through  IW31,  are  given  on  Table  1 and  are  shown  schematically 
on  Fig.  10.  Plenum  static  pressures  in  the  three  manifolds  are  denoted  PB01, 
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PB02,  and  FB03  for  the  outer  wall,  inner  wall,  and  base,  respectively.  Wall 
temperatures  were  measured  by  sixteen  thermocouples  located  on  the  inner  arid 
outer  walls  with  the  data  recorded  manually  on  the  thermocouple  readouts. 

The  locations  of  these  thermocouples  are  given  in  Table  1 and  are  shown  on 
Fig.  11.  These  thermocouples  are  labeled  TW01  through  TWl6.  In  addition, 
the  bulk  temperature  of  the  coolant  flow  was  measured  by  thermocouples  labeled 
TB01  through  TB06  as  shown  on  Table  1 and  Fig.  11. 

The  inlet  flow  distribution  was  measured  by  two  total  pressure  rakes  and 
two  wall  static  pressure  taps.  In  addition,  a measurement  of  circumferential 
distortion  was  made  through  four  midspan  total  pressure  rakes  located  at 
several  circumferential  stations. 

A midspan  total  temperature  probe  was  used  to  measure  the  bulk 
temperature  of  the  diffuser  flow.  The  location  of  all  the  inlet  flow 
instrumentation  is  shown  on  Fig.  12.  Diffuser  inlet  and  coolant  flow  rates 
were  measured  separately  by  orifices  located  as  shown  in  Fig.  12. 


Test  Results 

A summary  of  the  results  of  the  test  program  are  presented  in  the  test 
log,  Table  2,  and  a complete  set  of  test  results  in  Table  3.  Eleven  tests 
were  conducted  over  a period  of  two  days  at  three  simulated  power  settings 
and  three  different  coolant  flows  per  power  setting.  Two  tests,  2.01  and 
6.01,  are  considered  to  be  invalid  due  to  flow  perturbations  while  data  was 
being  recorded.  Each  test  was  repeated  to  provi  e valid  data.  The  main 
results  of  the  test  program  are  briefly  discussed  in  the  following  paragraphs. 

The  effects  of  film  cooling  on  pressure  recovery,  wall  pressures,  and 
wall  temperatures  were  evaluated  by  testing  coolant  flows  of  2.5,  5>  an<3 
percent  of  diffuser  inlet  flow.  The  largest  effect  of  this  coolant  flow 
variation  was  on  louver  wall  temperature,  as  shown  for  the  outer  wall  in 
Fig.  13.  Locations  of  the  outer  wall  cooling  louvers  are  shown  in  Fig.  7- 
Surface  temperatures  increased  rapidly  with  reduced  coolant  flow  and  the 
greatest  change  occurred  in  louver  4. 

The  effects  of  coolant  flow  on  wall  pressure  distributions  are  shown  in 
Fig.  14  for  a swirl  angle  of  zero  degrees.  Changes  in  wall  pressure  with 
coolant  flow  were  less  significant  at  swirl  angles  of  l6  and  21  deg  than  at 
zero  degrees.  Reducing  coolant  flow  rate  from  10  to  2.5  percent  on  inlet 
flow  resulted  in  a general  increase  in  wall  pressures  of  about  0.03  psi  for 
both  inner  and  outer  walls.  Note  that  the  increased  wall  pressures  were  not 
restricted  to  the  cooled  section  of  the  suppressor  but  also  extended 
upstream  to  the  inlet  static  pressures.  This  results  in  a decrease  in 
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pressure  recovery  (Cp)  from  0.63  to  O.58  with  coolant  flow  as  shown  in  Fig. 
15,  where  the  pressure  coefficient  Cp  is  given  by 


CP  = 


P-Pi 

P - P 
01  ri 


(139) 


However,  the  influence  of  coolant  flow  on  Cp  is  small  at  zero  swirl  and 
negligible  at  swirl  angles  of  16  and  21  deg.  Hence,  the  slight  ejector -effect 
of  the  film  coolant  at  zero  swirl  disappears  at  higher  swirl  angles. 


COMPARISON  OF  EXPERIMENT  AND  THEORY 


A set  of  calculations  were  made  to  demonstrate  the  capability  of  the 
computer  program  and  assess  the  code  by  comparing  the  theoretical  predictions 
with  experimental  data.  These  calculations  include  a solution  obtained  for 
turbulent  flow  through  a 6-deg  conical  diffuser,  termed  Fraser  Flow  "A" 

(Ref.  l8),  a baseline  calculation  for  the  ST9  demonstrator  diffuser  model 
tests  with  no  film  cooling,  and  a set  of  three  cases  for  the  flow  through 
the  ST9  demonstrator  IR  suppression  diffuser  corresponding  to  tests  logged 
in  Table  2.  The  calculations  were  numerically  stable;  however,  truncation 
errors  associated  with  the  linearization  of  the  differential  equations  were 
found  to  be  sensitive  to  streamwise  step  size.  The  sensitivity  to  step  size 
was  particularly  acute  for  the  slot  cooling  cases.  Since  the  slot  height 
was  only  one  percent  of  the  diffuser  inlet  height  in  these  cases,  the 
streamwise  step  size  had  to  be  very  small  (of  the  order  of  the  slot  height) 
to  allow  adequate  definition  of  the  resultant  flow  field  development.  Since 
the  computing  time  is  proportional  to  the  number  of  streamwise  stations  in 
the  calculation  field,  efforts  should  be  made  to  select  the  optimum  step  size 
in  order  to  minimize  computing  time  and  still  keep  truncation  errors  within 
reasonable  bounds. 


Fraser  Flow  "A" 


The  solution  for  the  turbulent  flow  through  a 6-deg  conical  diffuser  was 
obtained  with  the  current  calculation  procedure  and  compared  to  the 
experimental  data  of  Fraser  as  presented  in  Ref.  l8.  For  this  calculation, 
one  hundred  thirty  streamlines  and  fifty  streamwise  stations  were  used  to 
construct  the  coordinate  system  from  the  Schwartz -Christoff el  transformation 
as  shown  in  Fig.  l6.  The  mesh  distortion  parameter  placed  approximately 
twenty  mesh  points  between  0^Y+^10.  The  inlet  flow  was  constructed  from 
Coles'  velocity  profile  (Ref.  18)  with  the  inlet  boundary  layer  thickness 
taken  to  be  0.0528  in.  as  specified  by  the  Fraser  data  in  Ref.  l8. 


(l8)  Coles,  D.  E.,  and  E.  A.  Hirst:  Proceedings  Computation  of 

Turbulent  Boundary  Layers  - 1968  AFSOR-IFP-Stanford  Conference,  August  1968. 


As  demonstrated  in  Ref.  l8.  Coles'  profile  accurately  represents  the 
boundary  layer  mean  velocity  profiles  over  a wide  range  of  flow  conditions 
in  terms  of  two  parameters,  a friction  velocity,  U*  and  a wake  parameter, 
Coles'  profile  expresses  the  mean  velocity  distribution  by 

U+  z IT  ln  Y>  + 2 IT-  sin2  ( ) + B (i4o) 


where  U+  and  Y+  are  the  velocity  and  transverse  coordinate  written  in  wall 
variables 


U*  * U/U*  * U/Vt w //ow 


(141) 


(142) 

Y*  = YU%  « Y J T w/pw  /v 


tt  is  the  wake  parameter,  and  B is  a constant.  In  Eq.  (l4o)  the  first  term 
expresses  the  "law  of  the  wall"  and  the  second  term  expresses  the  "law  of  the 
wake".  Specification  of  U*  and  tt  uniquely  determines  the  velocity  profile. 

At  the  edge  of  the  boundary  layer  Y = 6 and  Eq.  (l4o)  becomes 


U 

U* 


2U 

K 


+ B 


(143) 


When  Eq.  (l40)  is  used  to  compute  the  boundary  layer  displacement  thickness 
6*,  the  relation 


S*Uco 

* 8u* 


= i + n 


(144) 


is  obtained.  Thus  given  5 and  6*,  the  wake  parameter  tt  and  friction  velocity 
U*  can  be  determined  from  Eqs.  (l43)  and  (l44)  and  then  the  velocity  distri- 
bution is  uniquely  determined  by  Eq.  (l4o). 
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It  should  be  noted  that  Fraser  Flow  "A"  is  a particularly  difficult  flow 
to  calculate  with  classical  boundary  layer  theory  which  assumes  a viscous 
flow  development  under  the  influence  of  a semi -infinite  nominally  inviscid 
outer  flow  field.  The  inability  of  classical  boundary  layer  theory  to  predict 
the  flow  field  development  is  shown  in  Ref.  18, where  a variety  of  boundary 
layer  theories  were  unable  to  predict  this  flow  field  accurately.  This 
inability  of  standard  boundary  layer  procedure  to  predict  Fraser  Flow  "A" 
stems  from  two  sources.  First,  the  wall  boundary  layer  reaches  the  diffuser 
centerline  approximately  halfway  downstream  from  the  diffuser  throat,  and  in 
this  region  no  potential  core  flow  exits.  Second,  the  Fraser  Flow  "A" 
diffuser  is  an  optimum  diffuser  in  that  it  is  designed  to  keep  a nearly 
separated  boundary  layer  as  the  flow  diffuses.  Boundary  layers  which  are  on 
the  verge  of  separation  are  very  difficult  to  predict  since  small  changes  in 
pressure  distribution  can  lead  to  large  changes  in  boundary  layer  thickness. 
However,  as  discussed  by  Coles  and  Hirst  (Ref.  18)  , the  Fraser  Flow  "A"  data 
are  reliable  and  based  upon  careful  development  of  a flow  configuration  with 
good  axial  symmetry.  Thus  the  flow  should  be  able  to  be  predicted  accurately 
by  a calculation  procedure  which  can  properly  account  for  the  disappearance 
of  the  central  potential  core.  The  present  analysis, which  does  not  assume 
the  existence  of  any  potential  core, yields  predictions  which  are  in  good 
agreement  with  the  data  as  shown  in  Figs.  17  and  18, where  comparisons 
between  experiment  and  theory  for  the  pressure  coefficient  and  wall  friction 
coefficient  are  present. 


Calculation  Procedure  For 
ST9  Demonstrator  IR  Suppression  Diffuser  Runs 


The  streamline  coordinate  system  used  to  predict  the  flow  in  the  FRDC 
diffuser  as  calculated  from  the  Schwartz -Christoff el  transformation  is  shown 
in  Fig.  19.  Since  this  diffuser  configuration  is  complicated,  it  is  helpful 
to  describe  it  in  some  detail.  First  it  is  noted  that  a short  "duct  inlet 
section"  has  been  added  to  the  diffuser  in  order  to  insure  that  the  initial 
inlet  flow  has  no  normal  pressure  gradient.  The  diffuser  centerbody  has  a 
blunt  "diffuser  base"  and  a lip  on  the  outerbody  (OD)  wall  extending  past  the 
base.  The  exit  centerbody  (ID)  and  OD  walls  were,  therefore,  extended  past 
the  diffuser  exit  plane  by  the  "extended  free  streamline"  show  in  Fig.  19. 

The  ID  and  OD  wall  contours  were  specified  at  fifteen  equally-spaced  stream- 
wise  stations.  The  computer  program  then  fitted  smooth  curves  through  these 
points  and  interpolated  the  curves  to  specify  eighty  streamwise  stations.  The 
input  mesh  points  are  indicated  by  a small  circle  and  the  diffuser  wall  by  a 
double  -ine.  In  addition,  seven  cooling  slots  are  located  on  the  wall  at  the 
stations  specified.  In  order  to  clearly  illustrate  the  location  of  the  slots, 
the  ID  and  OD  walls  were  opened  a small  amount  at  each  slot  so  that  at  the 
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exit  plane  the  wall  is  separated  from  the  wall  streamline  by  a small  amount. 
Finally,  six  struts  are  located  in  the  duct.  The  plan  view  of  one  strut  is 
indicated  in  Fig.  19  by  the  "strut  centerline",  "strut  leading  edge",  and 
"strut  trailing  edge". 

Eighty  streamwise  stations  and  one  hundred  thirty  streamlines  were 
calculated  for  the  streamline  coordinate  system.  The  mesh  distortion 
parameters  were  selectel  so  that  the  first  streamline  from  the  wall  was 
located  at  a distance  of  0.000077  times  the  duct  inlet  height.  This  placed 
approximately  five  mesh  points  between  0<Y+<10  at  a slot  exit  plane  and 
about  forty  mesh  points  within  the  slot  height . In  portions  of  the  flow 
removed  from  the  slots, approximately  thirty  mesh  points  were  placed  in  the 
viscous  sublayer.  For  all  cases,  the  inlet  flow  was  obtained  from  experi- 
mental data  shown  in  Fig.  20.  Instrumentation  in  the  inlet  plane  of  the 
diffuser  was  used  to  measure  total  pressure,  static  pressure,  swirl  angle, 
and  total  temperature  as  given  in  Table  3 for  each  test  case.  Twelve  data 
points  across  the  inlet  are  given.  These  data  were  then  interpolated  to  fit 
the  one  hundred  thirty  mesh  points  used  in  the  calculation, and  the  Mach 
number,  velocity  and  static  pressure, and  temperature  were  then  computed. 

This  experimental  data,  however,  is  not  accurate  enough  to  construct  the 
inlet  boundary  layer  profiles.  Therefore,  these  profiles  were  constructed 
from  Coles'  velocity  profile  by  specifying  the  displacement  thickness. 

Baseline  Case  - No  Film  Cooling  - No  Struts 

The  first  ST9  diffuser  case  calculated  consists  of  the  basic  diffuser 
with  no  struts  or  slot  cooling  and  is  termed  the  baseline  case;  the  baseline 
case  has  zero  swirl.  The  predicted  wall  static  pressure  distribution  is 
compared  with  the  experimental  data  taken  from  earlier  FRDC  cold -flow  model 
tests  (Ref.  7)  in  Fig.  21.  As  show.i  in  Fig.  21,  the  theoretical  predictions 
of  wall  static  pressure  distribution  without  cooling  and  without  struts  is 
higher  than  the  experimental  data.  This  ob  ervation  is  consistent  with  the 
observation  made  in  the  Fraser  Flow  "A".  In  addition,  it  is  noted  that  the 
analysis  predicts  the  appearance  of  a separation  bubble  on  the  centerbody 
wall  located  approximately  between  0.19  <Z/L  < 0.26;  this  prediction  is  in 
agreement  with  FRDC  model  test  data  which  also  shows  a separation  bubble  at 
the  same  approximate  location. 

A number  of  observations  can  be  made  about  this  test  case.  The 
predictions  of  the  analysis  are  in  qualitative  agreement  with  the  data  for 
pressure  coefficient  along  the  diffuser.  The  overall  pressure  rise  is  in 
fairly  good  agreement  with  the  data,  and  both  the  analysis  and  the  data  show 
a separation  bubble  in  approximately  the  same  duct  location.  However,  the 
detailed  variations  of  pressure  coefficient  is  not  in  good  qualitative 
agreement  with  the  data  and  in  the  region  of  the  separation  bubble,  the 
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analytical  predictions  and  the  data  show  significant  qualitative  discrep- 
ancies. However,  it  should  be  noted  that  it  is  very  difficult  to  maintain 
axisymmetric  flow  in  curved-wall  annular  diffusers  (see  Ref.  2).  This 
problem  is  compounded  when  separation  occurs  because  separated  flow  is  very 
sensitive  to  small  pressure  fluctuations.  Therefore,  the  size  and  location 
of  the  separation  may  be  affected  by  asymmetry  effects  which  m turn  may  lead 
to  discrepancies  between  prediction  and  data. 

STQ  Demonstrator  IR  Suppression  Diffusers  -2.5  Percent  Cooling  Rate 

The  2.5  percent  cooling  case  consists  of  the  ST9D  demonstrator  IR 
suppression  diffuser  run  with  six  struts  and  seven  cooling  slots  and 
corresponds  to  test  case  3-01  shown  on  Table  2.  This  test  case . represents 
the  diffuser  operating  at  60  percent  military  rated  power  and  film  cool.ng 
flow  rate  of  2.5  percent  of  the  diffuser  flow  rate.  The  coolant  plenum  total 
pressure  and  static  temperatures  are  given  in  Table  3.  Experimental  data, 
obtained  from  Table  3,  was  used  to  construct  the  inlet  flow  profiles  as 
described  for  the  Fraser  Flow  "A"  case. 

A comparison  of  the  predicted  flow  with  the  corresponding  experimental 
data  is  shown  in  Fig.  22  and  Fig.  23-  The  predicted  wall  static  pressure, 
distribution,  shown  in  Fig.  22,  indicates  a significant  drop  in  pressure  in 
the  region  of  the  strut.  This  effect  is  caused  by  the  strut  blockage  which 
accelerates  the  flow  through  the  strut  passage  (see  Ref.  7).  The  presence  o 
the  strut  exhibits  the  separated  region  which  is  present  in  the  baseline 
case,  and  separation  was  neither  predicted  by  the  analysis  nor  measured  in  the 
experimental  test.  The  first  slot  is  located  on  the  tip  wall  just  downstream 
of  the  strut  trailing  edge.  As  shown  on  Fig.  23,  the  wall  temperature  drops 
almost  by  700  deg  F at  this  slot  and  increases  rapidly  by  almost  500. deg 
as  rapid  mixing  of  the  film  cooling  flow  with  the  hot  diffuser  flow  is 
promoted . Successive  slots  follow  the  same  pattern.  The  predicted  wall 
static  temperature  distribution  is  in  good  agreement  with  experimental  data. 
The  wall  static  pressure  distribution  is  qualitatively  correct  for  ^/L<0.b, 
but  shows  significant  quantitative  discrepancies  with  the  data  Downstream 
of  the  fourth  slot,  the  wall  static  pressure  decreases  rapidly  due  to  large 
total  pressure  losses,  and  this  unrealistic  solution  is  not  plotted. 

STQ  Demonstrator  IR  Suppression  Diffuser  With  5 Percent  Cooling  Rate 

The  5 percent  cooling  rate  case  corresponds  to  test  case  2.02  shown  on 
Table  2.  This  case  represents  the  diffuser  operating  at  60  percent  military 
rated  power  with  a film  cooling  flow  rate  of  5 percent  of  diffuser  weight 
flow.  A comparison  of  the  predicted  results  with  experimental  data  is  shown 
on  Figs.  24  and  25.  The  wall  static  pressure  distribution  shows  very  li  - le 
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change  compared  to  test  case  3-01  with  2.5  percent  cooling  flow  rate  and  is 
consistent  with  earlier  observations  of  the  data  that  the  pressure  distribu- 
tion is  not  greatly  affected  by  the  coolant  flow  rate  as  shown  in  Fig.  14. 

The  wall  temperatures,  however,  show  a significant  reduction  from  a mean 
temperature  of  about  600  deg  F for  test  case  3*01  down  to  a mean  temperature 
of  about  450  deg  F for  this  test  case.  This  predicted  reduction  in  tempera- 
ture with  increasing  coolant  flow  rate  follows  the  trend  in  the  experimental 
data  shown  in  Fig.  13-  Again,  downstream  of  the  fourth  slot,  the  wall  static 
pressure  decreases  rapidly  due  to  large  total  pressure  losses  and  the 
solution  is  not  plotted. 

ST9  Demonstrator  IR  Suppression  Diffuser  With  10  Percent  Cooling,  Rate. 

The  10  percent  cooling  rate  case  corresponds  to  test  case  1.01  as  shown 
on  Table  2 and  represents  the  diffuser  operating  at  60  percent  military 
rated  power  with  a film  cooling  flow  rate  of  10  percent  of  diffuser  weight 
flow.  A comparison  of  the  predicted  results  with  experimental  data  is  shown 
on  Figs.  26  and  27.  Again,  the  wall  static  pressure  distribution  shows  very 
little  change  from  the  previous  cases,  indicating  that  the  coolant  flow  rate 
does  not  greatly  change  the  wall  static  pressure  distribution.  This 
observation  was  previously  noted  (see  Fig.  l4).  The  mean  wall  temperature 
shows  a drop  from  the  previous  cases  down  to  a mean  temperature  of  about 
300  deg  F and  demonstrates  the  expected  behavior  that  increasing  coolant  flow 
rates  reduce  wall  static  temperatures  as  shown  by  the  data  in  Fig.  13-  It 
should  be  noted  that  part  of  this  reduction  in  wall  temperature  comes  from  a 
reduced  plenum  total  temperature  as  seen  in  the  data  of  Table  3-  Therefore, 
an  accurate  calculation  of  an  IR  diffuser  performance  must  account  for  an 
overall  heat  balance  as  well  as  the  local  slot  cooling  effect. 


Discussion  of  Numerical  Calculations 

Since  this  report  presents  new  and  very  advanced  techniques  for 
calculating  turbulent  flow  in  ducts,  some  discussion  of  the  numerical 
problems  which  were  encountered  is  in  order.  In  particular,  suggestions  are 
made  on  means  to  improve  the  predictions  of  the  computer  code  for  slot  and 
film  cooled  problems.  The  Fraser  Flow  "A"  test  case  represents  a gooa  test 
case  for  the  purpose  of  checking  the  accuracy  and  reliability  of  the  computer 
program.  As  evaluated  by  Coles  and  Hirst  in  Ref.  18,  Fraser  Flow  A 
represents  reliable  measurements  based  on  careful  development  of  a flow  . 
configuration  with  good  axial  symmetry.  Pressure  distributions,  wall  fnc  ion 
coefficient  distributions,  displacement  thickness,  and  momentum  thickness 
distributions  are  presented.  In  addition,  accurate  measurements  of  the 
boundary  layer  profiles  are  given  for  eleven  axial  stations.  Finally,  he 
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Fraser  Flow  "A"  is  a clean  flow  showing  only  the  effects  of  boundary  layer 
growth  in  an  adverse  pressure  gradient.  Of  special  interest  in  this  est 
case  is  the  fact  that  the  flow  is  in  a nearly  separated  condition  for  a goo 
part  of  the  diffuser  length.  Complications  arising  from  struts,  slots,  swirl 

or  wall  curvature  are  not  present. 

Before  the  predictions  of  the  Fraser  Flow  "A”  case  were  compared  with 
the  experimental  data,  the  numerical  accuracy  for  the  computer  P-'°8™“ 
checked  for  internal  consistency  by  seveial  means.  The  most  important 
these  checks  was  a comparison  of  the  mean  flow  variables  obtained  by 
averaging  the  solution  for  each  dependent  variable  over  the  duct  height  and 
comparing  the  average  values  with  the  solution  for  these  same  variables 
obtained  by  integrating  the  one -dimen-tonal  mass  flaw  "lighted  ^erage 
euuations.  As  an  example,  these  equations  show  that  in  the  absence  o 
mLs  flow  bleed  and  wall  heat  transfer,  the  diffuser  weight  flow  and  mas. 
flow  weighted  average  total  temperature  are  constant.  An 

detailed  computer  printout  for  the  Fraser  Flow  A case  shows  that  ‘bleating 
variables  a-e  indeed  constant  to  at  least  five  decimal  places,  thus  indicat  np 
that  the  numerical  procedure  satisfies  the  integral  conservation  laws.  In 
refflrd  to  “e  prediction  of  Fraser  Flow  "A",  the  agreement  between  theory  and 
experiment  should  be  regarded  as  excellent  for  f in  motion  and  good  for 

the  pressure  coefficient,  particularly  in  view  of  the  fact  that  the  flow 
continuously  on  the  verge  of  separation.  As  shown  by  Stratford  (Be  f.  9). 
turbulent  boundary  layers  near  separation  may  take  on  a wide  _ 

profile  shapes.  The  shape  depends  upon  the  upstream  history  of  the  flow  a 
small^ changes  in  the  upstream  history  can  lead  to  large  changes  - -dary^ 
layer  development.  Thus  the  Fraser  Flow  A case  is  a diff  .. 
basic  analytical  procedure,  and  based  on  this  test,  1 is  cone,  u 
analysi^operates^well  in  predicting  flows  in  the  absence  of  such  complies- 
tions  as  struts,  swirl,  and  wall  curvature. 

Some  indication  of  local  errors  can  be  estimated  by  examining  the  total 
and  static  pressures  along  the  centerline  of  the  Fraser  Flow  A case.whio 
should  be  nearly  constant  since  the  viscous  and  heat  transfer  losses  are 
negligible  along  the  centerline.  This  comparison  indicates  local  accumulatea 
errors  of  the  order  of  0.1, which  is  quite  good.  From  these  resul  , - 

concluded  that  the  basic  numerical  procedure  is  accurate  and  that  any  improve 
ments  in  the  predictions  for  this  basic  case  must  come  from  a closer  exam 
tion  of  the  turbulent  mixing  length  model. 


(19)  Stratford,  B.  S.: 
Throughout  Its  Pressure  Rise. 


An  Experimental  Flow  With  Zero  Skin  Friction 
Journal  of  Fluid  Mechanics.,  Vol.  5,  1972, 


pp.  17-35 • 
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The  turbulent  mixing  length  model  may  be  an  important  source  of  error, 
producing  discrepancies  between  theoretical  predictions  and  experimental 
data.  Short  of  detailed  turbulence  measurements,  the  best  method  for 
assessing  this  effect  is  to  compare  predicted  boundary  layer  growth  with 
experimental  data  such  as  the  Fraser  Flow  "A"  test  case.  First,  it  should  be 
noted  that  the  Mach  number  for  the  Fraser  Flow  "A"  is  low,  so  that  0.1 
percent  error  in  static  pressure  represents  a 6.5  percent  error  in  static 
pressure  coefficient  (i.e.,  the  error  is  magnified  when  results  are  expressed 
in  terms  of  a pressure  coefficient).  Second,  it  is  noted  that  there  is  a very 
close  interaction  among  the  following  parameters:  mixing  length,  boundary 

layer  thickness  or  blockage,  and  static  pressure  gradient.  Thus,  increasing 
the  boundary  layer  thickness  decreases  the  static  pressure  because  of 
effective  blockage.  Increasing  the  pressure  gradient  increases  the  boundary 
layer  growth  because  of  the  additional  work  done  on  the  boundary  layer. 
Therefore,  a comparison  of  the  theoretically  predicted  pressure  distribution 
with  the  experimental  data  yields  an  indication  of  the  boundary  layer  growth 
and,  by  inference,  the  mixing  length.  Since  the  predicted  pressure  coefficient 
is  larger  than  the  experimental  data,  the  predicted  boundary  layer  thickness 
must  be  smaller  than  the  measured  boundary  layer  thickness.  Hence,  it  is 
concluded  that  the  mixing  length  is  too  small.  It  should  be  noted,  however, 
that  the  pressure  coefficient  prediction  is  very  good  up  to  a Z/L  = 0.4,  at 
which  point  the  boundary  layers  merge  (see  Fig.  17).  Furthermore,  downstream 
of  this  station  the  boundary  layer  is  in  a nearly  separated  condition  as 
shown  in  Fig.  l8.  Therefore,  it  is  concluded  that  modifications  of  the 
mixing  length  may  be  required  for  flows  with  merged  boundary  layers  or  nearly 
separated  boundary  layers. 

The  ST9  IR  suppression  diffuser  case  with  no  slots  or  struts  introduces 
an  additional  problem  in  making  an  accurate  flow  field  prediction.  This 
diffuser,  unlike  the  Fraser  Flow  "A"  diffuser,  has  significant  wall  curvature 
which  is  known  to  effect  the  mixing  length  (Ref.  20).  Thus  turbulence  can  be 
expected  to  increase  on  a concave  surface  and  decrease  on  a convex  surface, 
modifying  the  boundary  layer  growth  on  these  walls  and  changing  the  pressure 
distribution  accordingly.  Hence  it  is  expected  that  introducing  a better 
turbulence  model  which  accounts  for  wall  curvature,  would  produce  better 
predictions  for  wall  static  pressure  distribution  than  that  shown  in  Fig.  21. 
An  additional  complication  which  arises  in  the  baseline  calculation  is  the 
appearance  of  a separation  bubble.  Since  separation  is  a very  complex 


(20)  Bradshaw,  P. : Effects  of  Streamline  Curvature  on  Turbulent  Flow. 

AGAEDograph  No.  169,  1973. 
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phenomenon  which  is  very  sensitive  to  local  conditions,  predictions  of 
separated  regions  in  turbulent  flow  can  be  a strong  function  of  the. 
turbulence  model.  The  turbulence  model  used  in  the  present  effort  is  an 
equilibrium  model  based  on  measurements  of  turbulence  structure  in 
unseparated  flows  and,  therefore,  some  significant  error  may  be  introduced 
in  applying  a turbulence  model  based  upon  attached  turbulent  flow  data  to 
separated  flow.  Inaccuracies  of  the  turbulence  model  in  separated  regimes 
will  affect  the  predictions  downstream  of  flow  reattachment  since  total 
pressure  losses  produced  by  the  separation  bubble  cannot  be  reversed. 

For  slot  cooled  wall  cases,  test  case  3-01,  2.02,  and  1.01,  the 
comparisons  indicate  that  the  simple  mixing  length  model  used  in  this. report 
may  not  be  adequate.  The  turbulence  model  used  only  accounts  for  an  inner 
layer  mixing  length  influenced  by  the  wall  and  an  outer  layer  mixing  length 
influenced  by  the  free  stream.  The  mixing  layer  developing  between  the  cold 
slot  flow  and  the  hot  diffuser  flow  is  not  being  properly  modeled, and  this 
may  adversely  affect  the  theoretical  predictions.  If  the  turbulence  model 
used  is  not  appropriate  in  the  immediate  vicinity  of  the  slot,  then  any 
inaccuracies  would  be  compounded  for  diffusers  which  contain  a succession  of 
cooling  slots  such  as  the  ST9  IR  suppression  diffuser.  This  inadequate 
modeling  of  the  turbulence  structure  in  the  slot  mixing  region  may  explain 
the  increased  discrepancy  between  predicted  and  measured  pressure 
distribution  with  each  slot  (see  Figs.  22,  24,  and  26).  Thus,  in  summary, 
although  the  turbulence  model  appears  accurate  for  flows  in  the  absence  of 
curvature,  slot  cooling  and  separation,  it  may  need  to  be  refined  for  these 
effects  before  accurate  predictions  for  the  general  diffuser  flow  field  can 

be  made. 

The  numerical  method  used  in  this  report  can  be  subject  to  one  final 
important  source  of  error.  Since  the  linearization  of  the  equations  of  motion 
implies  a certain  degree  of  smoothness  in  the  solution,  any  local  errors 
introduced  in  the  initial  profile  may  produce  significant  errors  downstream. 
Usually,  as  in  the  diffuser  inlet  flow,  this  initial  error  is  dampened  and 
causes  little  difficulty  because  the  flow  variables  change  slowly.  At  the 
slot  interface,  however,  there  is  a large  temperature  discontinuity,  and 
since  the  flow  variables  (especially  temperature)  change  very  rapidly,  this 
initial  error  may  not  be  dampened.  Such  an  initial  error  would  cause  errors 
in  entropy  which  would  lead  to  inaccurate  predictions  of  the  pressure 
coefficient.  Indications  that  this  may  be  a factor  are  shown. by  the 
significantly  larger  errors  in  the  mass  flow  and  mass  flow  weighted  average 
total  temperature  which  are  not  as  well  behaved  downstream  of  a slot  as  they 
are  upstream  of  the  first  slot.  It  is,  therefore,  concluded. that  the  model 
used  to  predict  the  initial  flow  and  shear  for  each  slot  be  improved  so  as 
to  smooth  out  the  discontinuity  in  the  neighborhood  of  the  slot  exit  plane. 
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Based  on  the  experimental  data  presented  in  this  report,  it  is  conclude 
that  the  wall  temperatures  can  be  significantly  reduced  by  increasing  film 
coolant  rate.  For  coolant  flow  rates  less  than  or  equal  to  10  percent  of 
diffuser  flow  rates,  the  film  cooling  has  little  effect  on  pressure  distribu- 
tions or  pressure  recovery. 

An  examination  of  the  Fraser  Flow  "A”  test  case  demonstrates  that  the 
basic  numerical  methods  used  in  this  report  based  on  the  Schwartz -Christoffel 
transformation  to  calculate  an  orthogonal  coordinate  system  and  an  implici 
linearized  finite -difference  scheme  for  solving  the  equations  of  motion  for 
turbulent  flow  is  an  accurate  and  reliable  method  for  solving  internal  flows 
in  axisymmetric  ducts  of  arbitrary  wall  curvature.  Further  refinement  of  the 
turbulence  model  inregions  of  merged  boundary  layers  and  nearly  separated  flow 
is  indicated  by  the  comparisons  with  data. 

For  slot  cooled  walls  and  highly  curved  walls,  such  as  the  S19  IR 
suppression  diffuser,  further  refinement  of  the  turbulence  model  is  also 
indicated.  Specifically,  the  turbulence  model  should  include  the  effects  o 
wall  curvature  and  the  effects  of  a mixing  layer  between  the  hot  and  cold 
flows.  Finally,  the  initial  profiles  setting  up  the  slot  cooled  flow  need 
smoothing  of  the  temperature  and  density  discontinuity  in  order  to  minimize 
nonlinear  errors  in  the  calculation.  It  is  felt  that  if  the  indicated 
refinements  and  modifications  were  made,  the  resulting  computer  code  would 
have  a unique  capability  for  predicting  the  development  of  flow  fields  in 
axisymmetric  diffusers,  including  the  effects  of  wall  curvature,  struts, 
swirl,  and  film  cooling. 
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FIG.  3.  CONSTRUCTION  OF  SLOT  IN  DUCT. 


FIG.  6.  INSTALLATION  OF  SLOT  COOLING  SYSTEM. 
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FIG  8 SIDE  VIEW  OF  TEST  RIG. 


FIG.  9.  AFT  VIEW  OF  TEST  RIG 
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FIG.  12.  INLET  PLANE  INSTRUMENTATION. 
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PRESSURE  COEFFICIENT,  Cp  = (P  - Pi )/(Pq1  “ pl) 


FRASER  FLOW  "A"  CONICAL  DIFFUSER 


AFOSR  STANFORD  CONFERENCE  1968 
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FIG.  17.  COMPARISON  OF  EXPERIMENTAL  AND  PREDICTED  WALL 
STATIC  PRESSURE  DISTRIBUTION  FOR  FRASER  FLOW 
"A"  DIFFUSER. 
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WALL  FRICTION  COEFFICIENT,  Cf  = tw/(P01  -P^ 


FRASER  FLOW  "A"  CONICAL  DIFFUSER 
AFOSR  STANFORD  CONFERENCE  1968 


FIG.  18.  COMPARISON  OF  EXPERIMENTAL  AND  PREDICTED  WALL  FRICTION 
COEFFICIENT  DISTRIBUTION  FOR  FRASER  FLOW  "A"  DIFFUSER. 
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FIG.  19.  STREAMLINE  COORDINATES 
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ST9  DEMONSTRATOR  IR  SUPPRESSION  DIFFUSER  WITH  NO  STRUT 


DATA  TEST  NO.  3.01  0,60  MRP  SWIRL  ANGLE  0 DEG  '.OOLANT  FLOW  RATE  2.5% 


FIG.  20.  INLET  EXIT  MACH  NUMBER  DISTRIBUTION  . 
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FIG.  22.  COMPARISON  OF  EXPERIMENTAL  AND  PREDICTED  WALL  STATIC  PRESSURE 
DISTRIBUTION  FOR  ST9  IR  SUPPRESSION  DIFFUSER  WITH  2.5%  INJECTED 
COOLING  AIR. 


WALL  TEMPERATURE,  deg  F 
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TEST  NO.  3.01  0.60  MRP  SWIRL  ANGLE  0 DEG  COOLANT  FLOW  RATE  2.5% 
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FIG.  23.  COMPARISON  OF  EXPERIMENTAL  AND  PREDICTED  WALL 
TEMPERATURE  DISTRIBUTION  FOR  ST9  IR  SUPPRESSSION 
DIFFUSER  WITH  2.5%  INJECTED  COOLING  AIR . 
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ST9  DEMONSTRATOR  IR  SUPPRESSION  DIFFUSER 


TEST  NO.  2.02  0.60  MRP  SWIRL  ANGLE  0 DEG  COOLANT  FLOW  RATE  5.0% 


FIG.  24.  COMPARISON  OF  EXPERIMENTAL  AND  PREDICTED  WALL  STATIC  PRESSURE 
DISTRIBUTION  FOR  ST9  IR  SUPPRESSION  DIFFUSER  WITH  5.0°°  INJECTED 

COOLING  AIR. 
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FIG.  25.  COMPARISON  OF  EXPERIMENTAL  AND  PREDICTED  WALL  TEMPERATURE 

DISTRIBUTION  FOR  ST9  IR  SUPPRESSION  DIFFUSER  WITH  5%  INJECTED 
COOLING  AIR. 
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TEST  NO.  101  0.60  M-1P  CFR  = 0.10  SWIRL  ANGLE  0 DEG  COOLANT  FLOW  RATE  10% 


AXIAL  DISTANCE,  Z/L 

FIG.  26.  COMPARISON  OF  EXPERIMENTAL  AND  PREDICTED  WALL  STATIC  PRESSURE 
DISTRIBUTION  IR  SUPPRESSION  DIFFUSER  WITH  10%  INJECTED  COOLING  AIR 
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A* 
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Ak 
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bl 
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CD 


en0 

Ek 


Area  (ft^) 

Area,  a/r^  (dimensionless) 

Van  Driest  constant  (26.0) 

Critical  area  ratio  (dimensionless) 

Block  tridiagonal  matrix  (dimensionless) 

Diagonal  block  matrix  (dimensionless) 

Chord  (ft) 

Chord,  b/rr  (dimensionless) 

Location  of  pole  in  z plane  (dimensionless) 

Left  diagonal  block  matrix  (dimensionless) 

Speed  of  sound  (ft/sec) 

Drag  coefficient,  2D/(p2ufb)  (dimensionless) 

Friction  coefficient,  (dimensionless) 

Lift  coefficient,  2L/(p2U§b)  (dimensionless) 

2 2 

Specific  heat  pressure  ft  /(sec  deg  R) 

Pressure  coefficient,  (P  - \)/^01  - Pj_)  (dimensionless) 

2 2 

Specific  heat  volume  ft  /(sec  deg  R) 

Right  diagonal  block  matrix  (dimensionless) 

Drag/span  (lb/ft) 

Block  operator  matrix  (dimensionless) 

Streamwise  strain  (l/sec) 

Tangential  strain  (l/sec) 

Block  operator  matrix  (dimensionless) 
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®B 
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gb 

h+ 
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H+ 


LIJ 


m 


M 


0 

m 


m 


0 

M 


M 


Solution  matrix  (dimensionless) 

Force/area  (lb/ft  ) 

Force/span  (lb/ft) 

Complex  variable  source  solution  (dimensionless) 

Gap  between  walls  (ft) 

Gap  between  chord  lines  (ft) 

Gap  between  walls,  g/rr  (dimensionless) 

Gap  between  chord  lines  (dimensionless) 

Enthalpy  (ft2 /sec2) 

Height  of  inlet  duct  (dimensionless) 

Universal  stagnation  enthalpy,  (hQW"hQ)pwU*/qw  (dimensionless) 

Universal  adiabatic  stagnation  enthalpy,  (h0AW"hOA )/(u*)2  (dimensionless) 
2 2 

Entropy  ft  /(sec  deg  R) 

Element  of  L matrix 
Lift/span  (lb/ft) 

Matrix  for  k+1  point  (dimensionless) 

Mass  flow  (slugs/sec) 

Mass  flow,  m/(NBr^prUr)  (dimensionless) 

Mass  flow/area  slugs /(ft  sec) 

Universal  mass  flow  parameter,  R\^/(PwU  ) (dimensionless) 

Mass  flow/area,  m/(prUr)  slugs/(ft2  sec) 

Mach  number,  U/C  (dimensionless) 
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Nr 
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prt 
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*1 
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Q+ 

Qa 
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Ik 
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Rc 
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ft 


Streamwise  Mach  number  (dimensionless) 

Element  of  M matrix 

Boundary  condition  matrix  (dimensionless) 

Normal  coordinate  (dimensionless) 

Normal  coordinate,  n/(rrvr)  (dimensionless) 

Number  of  struts  (dimensionless) 

Reynolds  number,  rrprUr/4r  (dimensionless) 

Pressure  (lb/ft^) 

Universal  pressure  gradient  parameter,  p jy+'g  4x~  (dimensionless) 

Prandtl  number,  (-^  ) (dimensionless) 

Prandtl  number  turbulent,  ( )t (dimensionless ) 

Heat  flux,  -X  § 3b/ (ft  sec) 

Average  inlet  dynamic  pressure  (lb/ft2) 

Heat  flux,  q/(prUrCpTr)  (dimensionless) 

Universal  heat  flux,  q/qw  (dimensionless) 

Universal  heat  flux  (adiabatic),  q/(Pv^*^)  (dimensionless) 

Column  matrix  (dimensionless) 

Split  column  matrix  (dimensionless) 

Radius  (ft) 

Recovery  factor,  Eq.  (3.2.35)  (dimensionless) 

Radius,  r/rr  (dimensionless) 

2 2 

Gas  constant  ft  /(sec  deg  R) 
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r 
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St 
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T+ 

i^k 


'is 


un 
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UB 
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US 

Un 
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UB 


Radius  in  z plane  Eq.  (3*7.2)  (dimensionless) 
Radial  coordinate  (r,z)  plane  (dimensionless) 
Element  of  Rk  matrix 
Matrix  for  kth  point  (dimensionless) 

Strearawise  coordinate  (dimensionless) 

Strearawise  coordinate,  s/(rrVr)  (dimensionless) 
Stanton  number  (dimensionless) 

Blade  thickness  (ft) 

Temperature  (deg  R) 

Universal  temperature,  CpT/tf*2  (dimensionless) 
Column  matrix  for  k point  (dimensionless) 
StreauFise  velocity  (ft/sec) 

Normal  velocity  (ft/sec) 

Tangential  velocity  (ft/sec) 

Magnitude  of  velocity  (ft/sec) 

Blade  velocity  (ft /sec) 

Friction  velocity,  VTW/PW  (ft/sec) 

Strearawise  velocity,  Ug/Ur  (dimensionless) 
Normal  velocity,  Un/Ur  (dimensionless) 


Tangential  velocity,  U^/Ur  (dimensionless) 
Magnitude  of  velocity,  u/ur  (dimensionless) 
Blade  velocity,  UB/Ur  (dimensionless) 


u+ 

u* 

v 

V 

V- 

w 

W+ 

x 

X 

x 

y 

y 

Y 

X 


Y 


Y+ 


z 

Z 

Z 


Z 


B 


Z 


Universal  velocity,  U/U*  (dimensionless) 

Friction  velocity,  U*/Ur  (dimensionless) 

Metric  scale  coefficient  (dimensionless) 

Metric  scale  coefficient,  v/vp  (dimensionless) 

Volume  (ft3) 

Complex  variable  w plane  (dimensionless) 

Stream  function  inner  layer  (dimensionless) 

Distance  along  streamline  (ft) 

Distance  along  streamline,  x/rr  (dimensionless) 

Real  part  of  Z (dimensionless 
Imaginary  part  of  z (dimensionless) 

Distance  normal  to  wall  (ft) 

Distance  normal  to  wall,  y/rr  (dimensionless) 

Real  part  of  dw/dz  (dimensionless) 

Imaginary  part  of  dw/dz,  (dimensionless) 

Universal  distance  from  wall,  YpwU*/|iw  (dimensionless) 
Complex  variable  - z plane  (dimensionless) 

Axial  distance  (ft) 

Axial  distance,  z/rr  (dimensionless) 

Loss  coefficient  (dimensionless) 

Axial  coordinate  (r,z)  plane  (dimensionless) 

Column  operator  matrix  (dimensionless) 


1 

! 

! 
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Ens 

En0 


T1 

T1 
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e* 

0 
0* 

1 
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Swirl  angle  to  axis  (deg) 

Chord  angle  to  axis  (deg) 

Angle  Schwartz -Christoff el  transformation  (dimensionless) 
Ratio  of  specific  heats,  Cp/Cv  (dimensionless) 

Boundary  layer  thickness  (ft) 

Displacement  thickness  (ft) 

Boundary  layer  thickness,  6/rr  (dimensionless) 
Displacement  thickness,  6*/rr  (dimensionless) 

Streamwise  strain,  rrens/U r (dimensionless) 

Tangential  strain,  rren<lS/Ur  (dimensionless) 

Small  angle  in  z plane 

Transformed  normal  coordinate  (dimensionless) 

Imaginary  part  of  w (dimensionless) 

Blade/force  area,  rr,f/(pru^)  (dimensionless) 

Angle  of  streamline  to  axis  (deg) 

Momentum  thickness  (ft) 

Temperature  ratio,  T/Tj>  (dimensionless) 

Momentum  thickness,  Q*/lT  (dimensionless) 

Entropy,  (I-Ir)/R  (dimensionless) 

Von  Karman  constant  (0 .4l ) 

Thermal  conductivity  lb/(sec  deg  R) 


>1 


15  f 


fi  Viscosity  slugs/(ft  sec) 

H Plade  force/span,  ,f/(rrprU^)  (dimensionless) 

§ Real  part  of  w (dimensionless) 

n 3.14159 

II  Pressure  ratio,  p/pr  (dimensionless) 

p Density  (slugs/ft^) 

ps,pn  Radius  of  curvature  (ft) 

P Density  ratio,  p/pr  (dimensionless) 

Ps,Pn  Radius  of  curvature  (dimensionless) 

a Solidity,  b/gg  (dimensionless) 

£n,  Streamwise  stress,  Tns/(prU^)  (dimensionless) 

En ^ Tangential  stress,  Tn0/(prU^)  (dimensionless) 

Tns  Streamwise  stress  (lb/ft2) 

2 

Tn0  Tangential  stress  (lb/ft  ) 

t+  Stress,  t/t-w  (dimensionless) 

Column  matrix  for  boundary  conditions  (dimensionless) 
0 Tangential  coordinate  (radians) 

0C  Chamber  angle  (deg) 

0j  Angle  in  z plane  Eq.  (3.7*3)  (dimensionless) 

0B  Blade  dissipation  function  lb/(sec  ft  ) 

4>b  Blade  dissipation  function  (dimensionless) 
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X 

X 

0 

Y 


Clauser  constant  (0.0l6)  (dimensionless) 

Normal  coordinate  transform,  dT]/dn  (dimensionless) 
Stream  function  (slugs/ft) 

Stream  function  (dimensionless) 


Matrix  Operators 
T Transpose 


-1  Inverse 

Superscripts 

v Iteration  number 


Mean  or  average  quantity 

Variables  for  blade  force  calculation 

Deviation  from  mean  quantity 


Subscripts 

0 Stagnation  conditions 

Inlet  conditions 
Upstream  of  strut 
Downstream  of  strut 
Adiabatic 

Effective  turbulent 
Hub  conditions 
Incompressible  conditions 
Midspan  conditions 


1 

2 

3 

A 

E 

H 

I 

M 
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r Reference  conditions* 

T Tip  conditions 

W Wall  conditions 

•x'  Free  stream  or  edge  of  boundary  layer 


*Reference  conditions  are  based  on  standard  sea  level  atmosphere  conditions 
for  all  thermodynamic  quantities.  The  reference  radius,  rr,  is  the  inlet 
outer  radius,  and  the  velocity  is  the  mean  inlet  velocity. 
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